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We compute energy level correlations in weakly disordered metallic grains using the fermionic replica 
method. We use the standard cr-model approach and show that non-trivial saddle points, which 
break replica symmetry, must be included in the calculation to reproduce the oscillatory behavior of 
the correlations. We calculate the correlation functions in all three classical ensembles GOE, GUE 
and GSE including the finite-dimensional gradient corrections. Our results coincides with those 
obtained by the supersymmetric cr-model and the semi- classical trace formula. 
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I. INTRODUCTION 

The statistics of eneijgy levels of electrons in disordered metals has attracted much attention fiir-the past decades. 
Gor'kov and Eliashbergfcl conjectured that it obeys Wigner-Dyson-Jaws derived for random matricesH. This conjecture 
received a strong support almost twenty years later when EfetovtJ introduced the supersymmetric (SUSY) cr-model. 
It appears that the zero-dimensional version of the cr-model gives exactly the random matrix theory statistics of 
Wigner and Dyson. The use of the SUSY formulation seemed crucial since the alternative replica theoryoEl, when 
applied to the pure xandom matrix problem, seemed unable to reproduce the correct oscillatory behavior of the level 
correlation functiorJj. Later Altshuler and ShklovskiiQ realized that in a finite dimensional system, the correlation 
function is modified with respect to the universal random matrix level statistics: this modification appears when the 
energy difference is of the order of the Thouless energy, E c (equal to h over the diffusion time through the sample), and 
the corrections depend on the dimensionality d, conductance g, and shape of the sample. They used diagrammatic 
perturbation theory and consequently could trace only the modifications of the non-oscillatory part of the correlation 
functions. Finite dimensional modifications of the oscillatory part by the gradient terms were calculated in Ref.EJ, 
and more generally in RefB, using the SUSY technique. The result were subsequently reproduced using a semi- 
classical trace formula approach!!!]. It followed from these works that a power-law decay of the oscillatory correlations 
crosses over to an exponential decay at the scale E c . The precise behavior of this crossover depends both on the 
dimensionality and the symmetry class of problem. The essential feature of these results is that all non-universal 
terms may be expressed through the spectral determinant of a single classical differential operator. For the case of a 
disordered metal grain it turns out to be the diffusion operator in the corresponding geometry. 

In a recent paperELI we have shown how the fermionic replica method may be used to calculate the level statistics 
of the Gaussian unitary ensemble (GUE). The calculations of Refill were specific to the GUE and essentially used 
the Itzykson-Zuberlij integral for the unitary group. The purpose of the present paper is to present a more general 
approach to the fermionic replica calculations of the level statistics, which is not based on the peculiarities of the 
unitary ensemble, and uses rather the standard path of the cr-model. We shall present the calculations of the level 
correlations in disordered metals for all three classical symmetry ensembles: orthogonal (GOE), unitary and symplectic 
(GSE). We also include the effects of gradient terms on the level statistics reproducing exactly the results of Refs.ETEj. 
Our colleagues I.V. Yurkevich and I.V. Lerner have independently been developing a non linear cr-model approach 
with replica symmetry, using some complementary approach to oursE 3 ! n 

Our strategy is as follows: we deal with the standard fermionic replica cr-modelEl with an action written in terms of 
the (n\ + 112) x (n% + 712) dimensional Q matrix, where ni.2 — > are numbers of replicas. The symmetry group of the 
action G(ni +112) is broken down to the exact G(rix) x Gin^) by a finite energy difference u> = e% — 62 of the correlation 
function (G is a symmetry-group of—the Q-matrix, which depends on a symmetry class of the problem). Based on 
the experience of the GUE solutionis, we consider all possible saddle points of the cr-model both replica symmetric 
and replica non-symmetric. The latter spontaneously brake the exact symmetry of each sub-block G(n) down to 
G{p) x G(n — p) with < p < n (here n = n\, and p = p\,p2). The corresponding manifold of the Goldstone modes 
has an exact degeneracy for space-independent (q = 0) modes. The contribution of such saddle point manifolds (the 
coset space G(n)/G(p)G(n — p) ) to the partition (generating) function is proportional to their volume. The volumes 
of the coset spaces play a central role in our analysis, since after the analytical continuation n — > they determine 
which of the saddle points contribute to the generating function. It turns out that in addition to the replica-symmetric 
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(perturbative) saddle point (p = 0) there is only one additional saddle point manifold (in each block) with p = 1 in the 
GOE and GUE cases and two manifolds p = 1 and p — 2 in the GSE case. These replica-nonsymmetric saddle points 
give rise to an oscillatory part of the correlation function (since their action remains finite, and imaginary, in the 
n — > limit). One thus gets the correct oscillatory behaviors of the correlations, with one oscillation frequency in the 
GOE and GUE and two oscillation frequencies in the GSE. One should notice that the effect of the replica symmetry 
breaking saddle points are not limited to the correlation of levels. In the random matrix limit, they are known to 
describe the exponentially small tails of the density of states outside the asymptotic support of the spectrumli-3 and 
the finite size oscillatory correction insidelliJ. 

We then calculate the fluctuations around each of the saddle point manifolds in the Gaussian approximation. This 
is legitimate at relatively large energy w ^> A and for a good metal, g = E c / A 3> 1 (A is the mean level spacing 
and g is dimensionless conductance). No relation between u> and E c is assumed. As a result, one obtains the energy 
dependent amplitudes of the oscillatory parts, as well as those of the smooth parts, of the correlation functions, in 
the asymptotic regime ui 3> A. For small energy, lu <C E c , the correlation coincides with the random matrix theory 
predictions, whereas for largex-energy uj > E c it get modified in the non-universal (dimensionality and (/-dependent) 
way in agreement with Refs.Biil 

The structure of the paper is as follows. In section [h] we introduce notations and present a general discussion of 
the matrix field theory which allows to compute the energy level correlations. We compute the saddle points of this 
action and the quadratic fluctuations around them. We also comment on the connection of our approach to the usual 
non linear er-modcl. In section [II we apply the theory to the three classical ensembles. Finally in section [V we 
briefly discuss the results, their range of validity, and the possible further developments. The appendix contains the 
computations of the volumes of the relevant coset spaces, for each of the three symmetry classes. 

II. THE REPLICA MATRIX MODEL AND ITS SADDLE POINTS 

A. Preliminaries 

We shall discuss the correlation functions of the density of states (DOS), which is defined as 

v(e) = V- lr Tr8(ii + e-H), (1) 

where V is the volume and H = Hq + Udis is the Hamiltonian of the system. Here H$ is the Hamiltonian of the 
corresponding regular (clean) system, and Udis is a random disorder potential. We are interested in the large energy 
behavior and we thus measure all energies from the large positive chemical potential /i: the deviation e from [i is 
supposed to scale as the mean level spacing, A. 

The retarded/ advanced Green functions, G ± (e), are defined as 

G ± (e) = ( f i + e-H±ir ) )~ 1 (2) 

with rj - infinitesimal. The density of states v(e) is thus equal to the small r\ limit of [G~ (e)— G + (e))/(2-7ri). The average 
DOS at large enough p, is a featureless smooth function, which we shall approximate by a constant, (v{e)) = (Ay) -1 . 
Hereafter the angular brackets stand for the averaging over the ensemble of random disorder potentials, which we 
assume to be Gaussian and short-range correlated with zero mean and a variance given by 

(U dis (r)U di3 (r')) = (27ruT)- 1 S(r - r>) , (3) 

where r is an elastic scattering mean-free time. 

The main object of our study is the connected two-point correlation function of energy levels, defined as 

R{t 1 ,e 2 ) = V - 2 ( V {e 1 )v{e 2 ))-l. (4) 

Using the fact that (G ± G ± ) = (G ± )(G ± ) = -(-kVv) 2 , one finds 

i?(ei,e 2 )=2^2 (»5(ei,e 2 )-7r 2 ) ; S(e l7 e 2 ) = A 2 (G+( £l )G- (e 2 )> . (5) 



With the replica trick the two-point function S may be written ai 
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S(e u e 3 )= lim A ° * (Z^>^(E)), (6) 
m,2->0 n\Tii oeiO€2 
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where we have introduced the diagonal (m + n 2 ) x (m + n<i) matrix Eji — 5jiEj with 



E = ( /J, + e 1 +irj j = l,...,n 1 ; ,^ 

The generating function Z( ni > n2 > may be written as a functional integral over 2(ni + rWjJermionic fields. Getting from 
such a fermionic vector field theory to a matrix formulation is a standard procedure&El which we shall not repeat in 
details. Performing the Gaussian averaging over Udis, introducing a (ni +n 2 ) x (ni +n 2 ) Hubbard-Stratonovich matrix 
field Q(r) and integrating finally over the fermionic degrees of freedom, one obtains for the generating functionPEI 

Z^\E)= [d[Q}exp{-A[Q,E]], (8) 

where the action A[Q,E] is given by 

A[Q,E] = ^TrQ 2 -Tr^^E+^- + i-Q^ . (9) 
The symmetry of Q and the integration measure d[Q] depend on the symmetry class of the problem and will be 



discussed in section [III| separately for each ensemble. The trace operation includes both the replica indices and the 
spatial variables. 



B. Saddle points 

We shall evaluate the integral in Eq. (^) by the saddle point method, and check a-posteriori that such an evaluation 
is indeed justified in the limit of a weak disorder. A space independent solution of the saddle point equation satisfies: 

Qs. P . = — E^-T- 1 + f Q*-p) > ( 10 ) 

p x ' 

where the sum over p runs over the set of eigenmodes of the pure Hamiltonian, Hq, with some appropriate boundary 
conditions. The saddle point matrix Q s . p . may be diagonalised by a transformation: Q s . p . = U~ 1 AU where U is an 
element of the symmetry group G of the problem, U G G(n\ + n 2 ), and A is diagonal, Aji = SjiXj. The saddle point 
equation ([lO]) then implies that U takes the form: 

u =[ V o V 2 ) ' Vl G G( " l} ' V2 £ G( " 2) ' 
while each eigenvalue is a solution of the equation: 

i 1 

Xj = — V 5 : w sign (XA . (12) 

p ^1 ~ 2^ + 27 A J 

Here we have used the standard approximation, valid at large /it, where one substitutes the sum over modes by an 
integral over e p = p 2 /2m — /i, and neglects variations of v(t) in a vicinity of e = 0. We thus find that the eigenvalues 
of Qs.p. take the values: 

A J= ±1, j G l,...,ni +n 2 . (13) 

There exist (m + l){n 2 + 1) saddle point manifolds. Each such manifold M.(p lP2 ) may be indexed by two integers 
Pi G {0, ...,rti} and p 2 G {0, ...,n 2 }, and is generated by the block-diagonal symmetry transformations U of the type 
( pd| ) applied to the diagonal matrix: 

A (pip2) = diag{-l, ... - 1, +1, ... + 1, +1, ... + 1, -1, ... - 1} . (14) 

In order to scan the manifold M-( PlP2 ) in a non redundant way, one must restrict the symmetry transformations U to 
the coset space 
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G( ni ) G(n 2 ) 
"W) G(2 , l)G(ni _ pi) G(p 2 )G(n 2 P2 ) ■ U ° j 

It is useful to define the free propagator which is the diagonal matrix with eigenvalues 

V 2 * 



G i ^[E J + - + -X i ) . (16) 

The eigenvalues take four different values depending on the value of the index j. One can characterize them by two 
binary indices {a, a), where a = 1,2 designates two replica blocks with the energies e a and a = sign(Aj). For each 
energy (a — 1, 2), we have a retarded and an advanced propagator, G^, defined by: 

V 2 i x ~ 1 



G^^+- + £Q ±-j ; a=l,2. (17) 

In these notations the free propagator takes the form 

G = diag{Gr, ...G^,G+,...G+,G+,...G+,G 2 ,...G 2 }. (18) 



C. Saddle point action 

On the manifold M( PlP2 ) 1 the saddle point action is given by: 

V72 



A ( P1P2 ) = ^ Tr A 2 pip2) - Trln I E + j- + -A (pipa) ) , ( 10) 



where the trace involves both space and replica indices. In terms of the free propagators G a a defined in Eq. ( |l7| ) it 
reads: 

A (piP2) = ^Va"^ ~Pi TrlnG r - K "Pi)TrlnG+ -p 2 TrlnG+ - (n 2 - p 2 )Tr In G 2 , (20) 

where the traces involve only spatial variables. Expanding to the first order in re Q <C 1, and omitting unessential 
constant factors, one finds 

TrhiG± « T«7re Q /A (21) 
Finally neglecting all constants, which vanish in the limit ni 2 — > 0, one obtains for the saddle point action: 

A (piP2) = X ( ni£l ~ " 2£2 - 2pi£l + 2 P2£2) (22) 



D. Quadratic fluctuations 

Let us expand around the saddle point Q S . P . = U~ 1 A( piP2 - ) U, writing Q(r) — Q s . p . + J7 -1 <5Q(r)J7. The action 
expanded to the second order is diagonalised in terms of the Fourier components SQij (q) of the fluctuations: 



rii+n 2 

A[Q,S]«A (piP2) + - £ ^M ij (q)5Q ii (q)5Q ji (-q) (23) 



i,i=l q 

where the eigenvalue of the eigenmode (z, j, q), My(q), is given by: 



^) = ^A-i^E^(p+|)^(p-f) • (24) 
p 
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There exist a-priori sixteen different fluctuation eigenvalues for each momentum mode q. It is convenient to index them 
according to the binary decomposition introduced after Eq. (16). Each index j G 1, n + n' is associated with a pair 
of indices (a, a), where a = 1,2 characterizes the energy, e a , and a — sign(Aj) characterizes the retarded/advanced 
nature of the propagator. The sixteen different fluctuation eigenvalues are then: 



2rA At 2 



q 

P 2 



(25) 



The corresponding momentum sums are easily computed, resulting in: 

£<*(p+S)<*(p-§)-0i 



E^(p + ?)«5(p-f) 



2ttt 



I — Dq 2 r ± i(e a — e Q ')r 



(26) 



In the last expression we have expanded the sum to first order in the small parameters Dq 2 r <§; 1 and \ei — t2\T <C I 
where D is the diffusion constant defined as D = 2/xr/(md). 

We obtain eventually the following list of eigenvalues for each spatial mode q: 



When a = a' the eigenvalue is 



Mr. 



7T 1 

2A t 



We shall call the corresponding modes "massive" and denote their number for each spatial mode q as N„ 
• When a ^ a' and a/a' the eigenvalues are: 

7T 2 

M(l + )(2-) = M(2_)(l+) - ' n "~ 



(27) 



2A 

7T 



(Dq 2 - iw) 



M ( i_) (2 +) = M(2+)(i-) = ^{Dq 2 , 



'(2+)(i- 



2A 



(28) 



where to = e — e'. These are standard diffusive modes associated with the G(ni + n^) symmetry of the action, 
which is explicitly broken by a non zero u>. We shall call them "soft" modes and denote their number for each 
spatial mode q as Af s - and Af s + correspondingly. 



• When a ^ a' and a — a' the eigenvalue is: 



7T „ n 

2A 4 



(29) 



These are the Goldstone modes associated with the spontaneous breaking of the exact G(nx) x G(n 2 ) symmetry 
by replica non-symmetric saddle points. They exist only for the manifolds with non-zero p\ or We shall 
call the corresponding modes "zero" modes and denote their number for each spatial mode q as Af z ■ 

Such a separation of modes into massive, soft and zero is well justified in the limit where A <C oj; Dq 2 «l/r. This 
specifies the regime where our methods and results are applicable. The number of modes depends on the number of 
independent degrees of freedom of the Q matrix and should be specified separately for each of the ensembles. We can 
perform now the Gaussian integrals over SQ fluctuations. Each eigenmode with eigenvalue M^ arT ^ a i a i^ (q) contributes 
a factor 



Af( ao .)( Q ,/ (7 /j(q) 



(30) 



to the generating function (Z^™ 1 '" 2 '). The exception is the zero mode in the space independent, q = 0, sectorlij. 
This mode has identically zero mass, originating from the exact degeneracy of the M.( PlP2 ) saddle point manifold. 

Therefore, in the q = sector the integral over the zero mode results in the volume V^J^ of the coset space dl5|). 
These volumes are calculated in appendix for each of the three classical symmetry ensembles. 
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E. Generating function 



Finally putting all the factors together one finds for the average generating function 

AT s+ /2 / 2A N 



Pi ,P2=0 



n 



/ 2A 



V^q 2 



2A 



Dq 2 + iuj 



AT.+/2 



2A 



Dq 2 — iuj 



(2At)^/ 2 x 

M.-/2 



(2A-, 



\M m /2 



(31) 



where the saddle point action, ^4( PllP2 ) is given by Eq. (p2[). The first line in this expression represents the saddle point 
action and the fluctuations in the q = sector, whereas the second line originates from the Gaussian fluctuations 
of q ^ modes. In section [II we shall evaluate the coset space volumes, V/^^) an< ^ the number of modes, 
Nz,s,m = Nz,s,m(PiiP2)i for each of the classical symmetry classes. Hereafter we shall put A = 1, implying that 
all energies are measured in units of the mean level spacing, A. 



F. The non-linear sigma model 

Let us briefly comment on the connection to the usual formulation of the problem in terms of the non-linear a- 
model. For simplicity we discuss only the unitary case. This basically amounts to a reorganization of the computation 
we did above, which uses the strong hierarchy of masses (the massive modes are much more massive than the soft 
ones). Assuming first that 1/(tA) is large, one finds that the saddle points of Eq. (§1) are given by the set of matrices 
with Q 2 = 1. This set is actually an ensemble of m + n 2 + 1 disconnected manifolds S r , corresponding to all possible 
values of r = p\ — pi (or equivalcntly of the trace of Q). It is easily seen that all the modes which move away from 
these manifolds are massive, with a mass 7r/(2Ar). These massive modes correspond to perturbing the matrix Q by 
a SQ such that QSQ + 5QQ ^ 0, and the number of such massive modes is: n\ + n\ + 2r(r — m + n-i). Performing 
the integration over the massive modes one can write (up to irrelevant constants) : 

[rij+ri2+2r(r — n 1 +n 2 )]K 

Z^™\E)~ V (\/— I / d[Q]eM-A[Q,E}} , (32) 




s r 

where K is the number of different q modes. Notice that the manifold S r is characterized by Q 2 = 1, TrQ = 
n\ — n,2 — 2r. It thus contains all matrices of the type Q = U^ 1 A( pi P2 )U with U E G(n\ + n-i), and A.r pitP2 ) defined 
in Eq. (|4|), with pi — p 2 = T. 

Expanding the action for slow spatial variations of Q on the manifold S r , one gets to first order in u> the standard 
sigma model: 

[n\+nl+2r(r-ni+n2)]K 

IS. 



Z^\E)^ E V^T ' ^ / 5( i[Q]exp{-^f ( VQ) 2 -^Tr(AQ)}, (33) 



where A = A( .o) ■ 

For large to, one can study the sigma model by a saddle point approximation. The generic variations around a 
point Q of <S r , staying on 5 r , are of the type SQ = [Q, W], with an arbitrary matrix W(r). The stationarity of the 
action imposes Tr(A[Q, W)) = 0, which implies that the saddle points Q S . P . commute with A . One easily deduces 
that, on S r , the saddle points sub- manifolds are exactly the sub- manifolds M Pl . P2 with p\ — p2 = r. This approach 
basically reorganizes our previous computation by grouping together all the sub-manifolds with a fixed value of Tr Q 
(or pi — p2). As we shall see below, for ni,ri2 — ► 0, the only saddle points which contribute to leading order in At 
are the ones with pi = p 2 , which are all located on the same manifold with r = 0. Hence to the leading order one 
can approximate the generating function in (^) by an integral over the single manifold So, which is what is usually 
done in the cr-model approach. We control this result well at large u> because we can do the sums over the p\, pi 
and control the analytic continuation. But we believe that it is probably correct also for any cj. The reason is the 
following: in Eq. ( J33[ ) one may extend the sum over r to a sum going from — oo to oo, because when r is outside of the 
interval {—ri2, n\\ the volume of the integration space vanishes (this can be checked e.g. in the limit u> — * 0). One 
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may then take the limit n\, n 2 — > at fixed r. It is clear that the leading term comes from r = 0, which minimizes 
the number of massive modes. So the usual tr-model formulation, with an integral over Sq only and the action given 
by Eq. d33|), seems to be correct. However one must keep in mind that, on this manifold there are, for large ui, several 
saddle point sub-manifolds, which lead to the oscillations in the correlation functions. In the random matrix case, 
at least, (without the gradient term) one may also try to perform the integration over entire manifold Sq, without 
resorting to the saddle point method. If the non linear cr-model, formulated on So only, is indeed correct, this should 
give the exact result, not restricted tow> 1. This procedure was attempted in Rcf.El, but the analytical continuation 
of the expressions emerging from these calculations still remains to be studied. 



III. CORRELATION FUNCTIONS 



A. Unitary ensemble 



In the presence of a weak magnetic field the Hubbard-Stratonovich matrix Q is HermitianB. The corresponding 
symmetry group is the unitary group, G = U . The measure of the functional integral over hermitian matrices Qij(n) 
in Eq. (B|) is given by 



d[Q] = n 



II dQ " (q ) II (q) d5 ^ (q) 



1<J 



(34) 



There are (ni+n 2 ) 2 degrees of freedom for each spatial mode, q. Looking at the classification of modes, Eqs. ( p7j 
one finds that the number of massive, soft and zero modes is 

W m = p\ + (m - pi) 2 +pj + (n 2 - P2) 2 + 2pi(n 2 -p 2 ) + 2p 2 (n 1 -p\) = n\ +n% + 2(pi - p 2 )(pi ~Pi -ni + n 2 ) ; 
Af a+ = 2pip 2 ; 

W s _ = 2(ni -pi)(n 2 - p 2 ) ; 

Af z = 2p 1 (n 1 -pi)+2p 2 (ri2-p2), (35) 

which add up to (m + n. 2 ) 2 as they should. Notice that the number of zero modes, J\f z , coincides with the number of 
dimensions of the degenerate coset space manifold U {n\) / {U {p\)U {n\ — pi))x U (n 2 ) / (U (p 2 )U (n 2 — p 2 )) ■ The volume 
of this coset space is calculated in the appendix and is given by 



■>;(P1P2) _ (4 \Pi(n 1 -p 1 )+p 2 (n 2 ~p 2 ) pPi pP2 
k (nin 2 ) _ f «l r n 2 i 



where 



F p = 



r(i + n) 



r(i + P )T(i + n- P ) n r(l + (n - j + 1)) ' 



r(i + j) 



(36) 



(37) 



Since FP >n = the sums over p\ and p 2 in Eq. ( |3l|) may be extended up to infinity. The resulting expression may 
be then continued analytically to n\. 2 — > 0, using the procedure which was detailed in our previous paperO. The 
analytical continuation, n — > at fixed p, of the F% symbol, Eq. (pTu, shows that 



F° 



1: F 1 



n- F p - 2 



0(n p ) . 



(38) 



Therefore only the terms with pi )2 = 0, 1 may contribute to the correlation function S, Eq. (^). The number of 
massive modes, in the limit nx,i — > at fixed pi,p 2 , is J\f m —> 2(pi — p 2 ) 2 . Therefore the terms with p\ 7^ p 2 can be 
neglected to leading order in the parameter At <C 1. One thus ends up with the two contribution to the generating 
function: p\ = p 2 = and pi = p 2 = 1. 

The first piece with p\ = p 2 = is the usual replica symmetric contribution. Using Eqs. (|3T|), ( |35| ) and ( p8[ ) one 
finds: 



Dq 2 



ito 



(39) 



Using Eq. 



one finds for the corresponding contribution to the correlation function: 
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(40) 



This is the well known perturbative contribution!]. The replica non-symmetric manifold with p\ = p2 = 1 gives: 

2ttiw / (£)q 2 ) 2 



q#0 



(Dq 2 ) 2 +uj 2 



(41) 



Differentiating over ei and e 2 according to Eq. (^|) and keeping only the leading contribution in uj/A 3> 1, one obtains 
for the corresponding contribution to the correlation function 



S(w)| Pl = P2 =i - - 

where we have introduced the spectral determinant of the diffusion operator, defined as 



z>m = n 

q/o L 



Dq 2 



Finally using Eq. (H) one finds 



„/ \ 1 „ x—v 1 cos(27toj) „ 



(L>q 2 - iw) 2 2tt 2 w 2 



(42) 



(43) 



(44) 



in agreement with Refs.BII 



B. Orthogonal ensemble 

If the time reversal symmetry is not broken the Hubbard-Stratonovich matrix Q appears to be a self-dual real- 
quaternion matrixd. This means that each element may be written as 



3 

E< 

a=0 



with real Qf, , where 



to = 



1 
1 



T\ = 



T2 = 



T3 



-i 
i 



(45) 



(46) 



Moreover Qji = (Qij) •> where conjugation operation acts on the r matrices. Such matrices may be diagolized by 
rotations from the symplectic groups, G = Sp{n\ + na), which is the relevant symmetry group for the GOE. Diagonal 
elements of the Q-matrix, Qu are characterized by a single real number, whereas off-diagonal ones Qi<j are 
parameterized by four numbers, Q^, a = 0, . . . , 3. Altogether there are 2{n\ + n 2 ) 2 — {n\ + n 2 ) degrees of freedom 
for each spatial mode, q. The measure of the functional integral in Eq. (pi) is given by 



= n 



IK^nil^tq) 

j i<j a=0 



(47) 



One easily finds that the number of massive, soft and zero modes is 



(ni - 
4(pi 



M m = 2p{-p 1 + 2(n 1 ~ Pl )'- 
= 2n\ + 2«2 — ( n i + n 2) 
M a+ = 4pip 2 ; 

M 8 - = 4(m -pi)(n 2 -P2) ; 

A4 = 4pi(ni - pi) + 4p 2 (n 2 - Pa) , 



pi)+2p|-p 2 + 2(n 2 
-p 2 )(pi -pa-ni +n 2 ) 



P 2 ) 2 - (n 2 -P2) +4pi(n 2 -p 2 ) + 4p 2 (ni - pi) 



(48) 
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which add up to 2(ni + t^) 2 — (ni + 712)- The number of zero modes, Af z , coincides with the number of dimensions of 
the degenerate coset space manifold Sp(ni) / (Sp(pi)Sp(ni — pi)) x S , p(n 2 ) / (Sp(p2)Sp(ri2 — P2)) ■ The volume of this 
coset space is calculated in the appendix and is given by 

y(PiP2) _ (4 \2pi(ni-pi)+2p 2 («2-P2) ftPl £P2 (49) 
(nina) ^ ' n i n 2 ' V ' 

where 

r(i + n) A r(i + 2j) 



n m+ (50 » 



r(l +p)T(l + n-p) 1 = 1 r(l + 2(n - j + 1)) 



Since -E p> " = the sums over p\ and P2 in Eq. ( 
then continued analytically to 711,2 (cf. Ref. 



may be extended up to infinity. The resulting expression may be 
. In the limit n — > the symbol, Eq. (50), is given by 



£?Uo = l; E^ = 2n; = 0(n*) . (51) 

Therefore only the terms with p\_2 = 0, 1 may contribute to the correlation function S(u>), Eq. (^). The number of 
massive modes, in the limit where 711,2 — *■ at fixed pi,p2, is Af m — * 4(pi — j>2) 2 - Therefore the terms with p\ 7^ P2 
may be neglected to leading order in the parameter At <C 1. Like in the unitary case only two terms with p\ = P2 = 
and pi = P2 = 1 contribute to the generating function. 

The replica symmetric contribution p\ = p 2 = is very similar to the one of the unitary ensemble. Using Eqs. (|3l|), 
(fl8|) and d|l]) one finds: 

(, \ 2niri2 

^— j ■ (52) 

Using Eq. (Q), one finds for the corresponding contribution to the correlation function: 

%)U =0 ^ 2 +2^ (i)q2 1 ^ , (53) 

in agreement with the known perturbative calculations!. The replica non-symmetric saddle point manifold, p\ = 
P2 = 1, contribute as 

Differentiating over t\ and £2 according to Eq. (^]) and keeping only the leading contribution in uj/A 3> 1, one obtains 
for the corresponding term in the correlation function 

p 2iriuj 

SH| Pl=P2= i = ^2? 2 M, (55) 

where the spectral determinant, T>(u>), is defined by Eq. (|43|). Finally, from Eq. (|^) one finds 

1 ~, 1 cos(27tcj) o , . . 

*M = ^E w _^ )2 + ' (56) 



again in agreement with RefsBEEl 



C. Symplectic ensemble 

If the spin of electrons is taken into account and the strong spin-orbit scattering is assumed the Hamiltonian of the 
system acquires a quaternion (symplectic) structured. The corresponding symmetry of the Q-matrix is the orthogonal 
one, G = 0(ni + n 2 ). The Q is a real symmetric matrix and the integration measure in Eq. (pi) is 
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d[Q\ = n 



i<j 



(57) 



There are [(ni + n 2 ) 2 + («i + n 2 )]/2 real degrees of freedom for each spatial mode, q. The number of massive, soft 
and zero modes is 

M m = ^[pf +Pi\ + -Pi) 2 + ( n i -Pi)] + ^[p 2 + Pa] + \l( n z -Pz) 2 + ( n 2 -Pa)] + Pi(n 2 -pa) +Pa(ni -pi) 

= ~[nf + nl + n! + n 2 ] + (pi - p 2 )(pi -p 2 — n\ + n 2 ) ; 
■M>+ = pip 2 ; 

TVs- = (m -pi)(n 2 -p 2 ) ; 
■A4 = Pi(«i - Pi) +Pa(«a - Pa) , 



(58) 



which correctly add up to [{n\ + n 2 ) 2 + (n± + n 2 )]/2. The number of zero modes, Af z , coincides with the number of 
dimensions of the degenerate coset space manifold 0(ni) / (O(pi)O(ni — pi)) x 0(n 2 ) / ' {0(p 2 )0(n 2 — p 2 )) . The volume 
of this coset space is calculated in the appendix and is given by 



y(PlP2) _ (2 /p)Pl(™l-Pl)+P2("2-P2) QPl QP2 

(niri2) v v / m n 2 ' 



(59) 



where 



Gl 



r(l + n) 



r(i+i/2) 



r(i + P )r(i + n - p) r(i + (n - j + 1)/2) 



(60) 



Since G^ >n = the sums over pi and p 2 i n rEfl. ( |3l| ) may be extended up to infinity. The resulting expression may be 
continued analytically to n ii2 — * (cf. Ref.liil). The G p n symbol, Eq. (Eft]), in the limit n — > is 



-n ; G n ^ 



-ft, <-r„^o 



0(n 



I(p+1)/2] 



(61) 



where [a;] denotes integer part of x. Therefore only the terms with pi j2 = 0, 1,2 contribute to the correlation function 
S{u)), Eq. (§). The number of massive modes, in the limit where ni :2 — + at fixed pi,p 2 , is M m — > (pi — p 2 ) 2 , making- 
terms with pi ^ p 2 small in the parameter At <C 1 . One therefore finds three relevant contributions to the generating 
function: p x — p 2 — 0, p\ — p 2 = 1 and Pi = p 2 = 2. 

The replica symmetric contribution pi = p 2 = comes almost without changes. Employing Eqs. ([5l]), ( |58| ) and 
( |6l| ) one finds: 



1 



£>q 2 — ioj 



ni« 2 /2 



(62) 



From Eq. (g), one finds for the corresponding contribution to the correlation function: 

1 



S(uj)\ Pl =p 2 =0 = 7T 2 + 



(Dq 2 - ilj) 2 



(63) 



in agreement with Ref.Ll. The first replica non-symmetric manifold, p% =p 2 = 1, results in 



in 2 ■ 



8a; 



n 

q#0 



( (Dq 2 ) 2 



V(Z?q2 



1/2 



(64) 



Differentiating over ei and e 2 according to Eq. (^|) and keeping only the leading contribution in ui/A 3> 1, one obtains 
for the corresponding contribution to the correlation function 



S(u>)\ Pl 



Pl=P2 = l 



2w 



(65) 
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where the spectral determinant, T>(uS), is defined by Eq. (43). Finally, the second replica non-symmetric manifold, 
Pi = P2 = 2, gives 



=2 = nin 2 



n 



( W) 



2\2 



(4tt) 4 u; 4 11 V(£»q 2 ) 2 + lo 2 



and consequently 



Using Eq. (0) one finds 



SHI 



Pl=P2=2 



47ritJ 



16tt 2 w' 



cos(27Ttj) , — — — - cos(47rw) 



47r 2 ^ (Dq 2 - iu) 2 Abj 



32ir 4 uj 4 



(66) 



(67) 



(68) 



again in agreement with Refslrc3 



IV. DISCUSSION OF THE RESULTS 



Let us briefly discuss the energy scales, the approximations involved in the calculations and their range of validity. 
There are four important energy scales: the mean level spacing, A; the Thouless energy, E c = %D/L 2 (L is the 
system size); the inverse scattering time, tt/r; and the chemical potential, fi. In the calculations above, the following 
hierarchy was assumed: A <C E c <C h/r <C /i. The condition h/r <C M, which means that the disorder is weak 
enough, was used to evaluate momentum sums by contour integration. The inequality E c <C tt/r, which is equivalent 
to L much larger than the mean free path I, tells that the system is in the diffusive regime. It was used to derive 
the diffusive dispersion law in Eq. (|2^). Finally, g = E c /A is the dimensionless conductance, and the condition that 
j> 1 means that the system is metallic. This condition was used to calculate integrals over zero modes with 
in the saddle point approximation. One more inequality was assumed in our derivation, the fact that the difference 
in energies uj is much larger than the level spacing A. This is a technical assumption, which allowed us to evaluate 
soft modes integrals by the saddle point technique. It would be interesting to perform the calculations without this 
last assumption, extending thus the results to arbitrarily small lu. 

Our calculations give the correlation as functions of u) in the form of a finite sum of oscillating harmonics (two in 
the GOE and GUE and three in the GSE), with u> dependent amplitudes. The set of harmonics is exact and has to do 
only with the symmetry of the problem, specifically with the volumes of the relevant coset spaces. The amplitudes, 
on another hand, were obtained in the saddle point approximation only. Using our formulation, one may develop a 
perturbation theory near the replicsupon-symmetric saddle points, much in the same way as it was done near the 
replica symmetric one, see e.g. Ref.tS. From such a perturbation theory one may obtain a systematic expansion of 
the amplitudes of the oscillatory terms, in powers of A/ui < 1. n 

We would like to point out striking similarities between our replica approach and the SUSY one of Ref E which was 
also based on the saddle point calculations. In particular the list of modes is the same. In the SUSY case, the zero 
modes and soft modes are respectively associated with the rotations inside the fermionic block and between fermions 
and bosons. In some sense our p and n — p replica blocks are similar to the bosonic and fermionic blocks of the 
SUSY theory. To appreciate better this analogy, one would need a more detailed understanding of the mathematical 
structure of the theory. In particularly, one would like to define the unitary (or other) group, U(n), for non-integer 
n and trace its relation to the graded symmetry. Another interesting problem is to appreciate better connections to 
the semi-classical method of Ref.H3. 

The existence of the replica non-symmetric saddle points opens two very important questions. One concerns their 
relevance to the renormaiization group treatment of the localization pxoblem for one electron. Another, even more 
challenging one, is to extend the replica theory of interacting electronstZl to account for new saddle points. 
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APPENDIX: ZERO MODES AND VOLUME OF THE COSET SPACE 



A. Unitary case 



The manifold M( PlP2 ) of saddle point matrices Q is generated by unitary transformations U of U(ni) x [/(712), 

applied to the diagonal matrix A( piP2 \, cf. Eq. We must first find which choices of V\, V2 actually change the Q 

matrix. A general unitary transformation V of U\n-\) can be written as a product RW where the matrix W has the 
block diagonal structure: 



W = 



W P1 

W n 







(69) 



where W Pl and W ril - Pl are unitary matrices of size p\ and n\ — p\ respectively. The matrix W belongs to the 
subgroup U{pi) x U(ni —pi) of U(n\) which leaves the saddle point matrix invariant. The 'proper' V\ transformations 
which change the matrix Q while staying on the saddle point manifolds are thus the elements R\ of the coset space 
u( P ^u(n 1 - Pl ) ' ano - similarly the proper V2 transformations are elements R2 which belong to u( P ^uln 2 - P 2) ' 

To compute the volume of the sgt of proper R transformations in U(n) (here n stands for either n\ or 112), we 
start from the usual decompositions of the integral over the group ofnxn Hermitian matrixes X in terms of the n 
eigenvalues Xj, and the unitary transformation V such that X = V^ _1 (diag{xi, ...x n })V: 



dp n (X) = dxj Y[ 



x i) n 

l<j<k<n 



x k ) 2 dp n {V) 



(70) 



In this integral we have ordered the eigenvalues (the 9 function is Heavyside's step function), in such a way that the 
integral over V scans the whole set of allowed unitary transformations. We can compute the normalization of the 
'angular' measure for instance by integrating a Gaussian function: 



dp n (X) exp 



-TrX : 
2 



_ 7r n 2 /2 2 «/2 



(71) 



which can be computed using the above decomposition and the Selberg's integral! 



n! 



/ d Pn (V) J f[ dx, I] {xj - x k f exp - 1 x n = ^ / d Pn W 



3=1 j<k 

Therefore one gets the normalization of the integral over the angular measure: 

n ! 



(2^ /2 I]r(j + i) 

3=1 



V 



u 



n-=ir(j + i) 



(72) 



(73) 



This result is easily checked by a direct counting argument: the choice of V is a choice of a Hermitian basis. The 
first vector of the basis is an arbitrary unit vector, the corresponding volume of integration is thus 6*2,1/271 where 
Sd = dit d / 2 /r(l + d/2) is the volume of the d dimensional unit sphere, and the division by 2n deals with a global 
phase choice. The second unit vector of the basis must be orthogonal to the first one, which fixes two real conditions, 
and its volume is thus Sm-il^-K. After iterating this construction, one gets the result (f73|). 
We now decompose V = RW, and the angular integral dp n (V) as: 



dp„(V) = dp p (W p )dp n - p (W n - p )dp ntP (R) 



(74) 



This defines the measure dp n ^ p in the 2p(n — p) space of the proper transformations R. The normalization of this 
measure is: 



J dp n ,p(R) = 

where we have introduced the symbol F% defined by 



y II 



V U V U 



(75) 
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We can now go on to the exact evaluation of the zero mode integrals. We keep within the subspace of the space 
independent Q matrices (q — modes) which are the only modes having a zero eigenvalue sector. Clearly the zero 
modes integrals factorize into two independent pieces, associated with each of the two coset spaces ;7(pi)t7(ni-pi) an d 



Cg,lCL10 1CXI_. IU1 IZjVJ _VVLJ llll_Hjpi_.llV.HJl.LL J_»HJL_,(Ji_, CLOOVJ^ICL UGL-l VVLL11 CCL^ll VJ1 LL1C LWU *~\J OLJL DJJCL,L_UO jj JT 

^(" 2 ) We can compute each such piece by working with a n x n Hermitian matrix X and computing: 



U(p 2 )U(n 2 -p2) ' 

= J d Pn (X) exp (-^TrX 2 + Trln (e + ^-X^j (77) 

We expand around the saddle point manifold generated by X = A p = diag{Ai...A„} = diag{— 1..., — 1, 1, ...1} by 
writing: 

X = RApR- 1 +SX , 8X = W(diag{xi, Xn^W^RT 1 , (78) 

where W, as above, is in U(p) x U(n — p) and J? is a proper transformation. Using the decompositions of the measure 
defined in Eqs. (|70|), (ff4|), one obtains 



Z {n) = I dx x ...dx p JJ { x i~ x if Yl ( x j+i- x j) [dp p (Vp)x 

l<i<i<p 3=1 J 

/n— 1 /. 
dx p+1 ...dx n JJ (xi-Xj) 2 Yl 0(x J+ i~x j ) I dp n - p (V n -p) x 

p+l<i<j<n j — p+1 

p P n 

dp rhp {R)Yl Yl {-2 + x t -x j ) 2 0(2 + x p+1 -Xp) x 
J i=ij= p +i 

exp S^-A p + ^ J2 x ) + ^2 Tr H E + «A p /2r)- 1 «(_B + _A p /2r)- 1 <5X) | 



(79) 



where A p is the saddle point action: 

A p = — TrA 2 + Tr In (E + iA p /2t) . (80) 
4r 



The integral in Eq. ( |79[) can be simplified by the following observations: the integrals over xi are all massive modes, 
and thus one can assume that \xi — Xj\ <C 1. Therefore the third line of Eq. ([79]) is just a constant, equal to 
7r p(n-p)^7'P22p(n-j)) Apart from this constant, the rest of Eq. ( f79| ) is nothing but the integrals over the massive modes. 

What we have shown here is that, in the sector q = of uniform fluctuations, the exact integral over the saddle 
point manifold (the zero mode directions) gives a factor: 

( 47r )Pl(«l-Pl)+P2("2- P 2) pPl pP2 ( gl ) 



B. Orthogonal case 



We shall not repeat here all the steps of the previous computation, they run in exactly the same way. We just give 
the main modifications. The integral over the symplectic group, generalizing Eq. ( |72] ) is equal to: 



V* = / dp n (V)=n\ 



d Pn (X) exp (-^jij JJ dXj JJ ( 



xj — Xk) exp 



j 



(82) 



and the computation of Selberg's integral gives the volume: 
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The ratio of volumes is 

V s 

y n 



"p "n—p 



where 



™ = nl rg =1 r(2j + l) ng£(2i + l) 

^- P !(n-p)i rrur(2 J + i) • 

Finally the factor coupling the eigenvalues j < p to those j > p in the analog of Eq. ((79|) becomes: 

p n 

n n {-i+xi-xjf (86) 

so that the final integral over the zero mode manifold is: 

^4 7r ^2[pi(ni-pi)+p 2 («2-P2)]£jPi £P2 ^87^ 



C. Symplectic case 

We just give again the main modifications. The integral over the orthogonal group, generalizing Eq. (|72[) is equal 



to: 



V° = / d Pn (V) = n! / dp„(X) exp 



2 



/ n n - ^1 ex p 



j=i j<k 

where X is a real symmetric matrix. The computation of Selberg's integral gives the volume: 

7r n 2 /4+n/4)2« 2 /4-5n/4 



V?=n\- 



The ratio of volumes is 



where 



r n 



V°V° 



m =1 r(i+j/2) 



( 27r )P(»-P)/2 G P 



Gf, = 



i\ n?=ir(j72 + i) n;d p r(j72 + i) 



p!(n-p)! n"=ir(i/2 + l) 

Finally the factor coupling the eigenvalues j < p to those j > p in the analog of Eq. (^9h becomes: 



p n 



n n i -2+^-^-1 

i— 1 j— p+1 

so that the final integral over the zero mode manifold is: 

(4 7r )[Pl(ni-Pl)+P2("2-P2)]/2 QPl QP2 



(89) 



(90) 



(91) 



(92) 



(93) 
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